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Abstract 

The charge radii of octet-baryons obtained in quenched lattice-QCD calculations are extrapo- 
lated within heavy-baryon chiral perturbation theory. Finite-range regularization (FRR) is applied 
to improve the convergence of the chiral expansion and to provide estimates of quenching artifacts. 
Lattice values of quark distribution radii and baryon charge radii for in the range (0.1, 0.7) GeV^ 
are described very well with FRR. Upon estimating corrections for both finite-volume and quench- 
ing effects, the obtained charge radii of the proton, neutron and T,~ are in good agreement with 
experimental measurements. The predicted charge radii of the remaining octet-baryons have not 
yet been measured and present a challenge to future experiments. 
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I. INTRODUCTION 



The study of the electromagnetic form factors of baryons is of crucial importance to un- 
derstanding the non-perturbative properties of QCD. Though QCD is accepted as the funda- 
mental theory of the strong interaction, it remains a theoretical challenge to quantitatively 
probe the non-perturbative domain. There are many effective methods and phenomenologi- 
cal models which have been applied to study the electromagnetic properties of baryons: the 
cloudy bag model [H, the constituent quark model [1,01, the l/N^ expansion approach 4|, 
the perturbative chiral quark model |5|, the extended vector meson dominance model 6|, 
the quark-diquark model 0] and the Schwinger-Dyson formalism [sl, 0, 10 1. Various formula- 
tions of heavy-baryon chiral perturbation theory (xPT) have also been widely applied to this 
problem 0, El, El, El . It has been observed that expansions in xPT are consistent with 



experimental results up to Q — 0.1 GeV [Ij]. Extensions of xPT to explicitly incorporate 
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vector mesons have been demonstrated to improve the applicability to ~ 0.4 GeV^ 

As well as the above model calculations, the past few years have seen increased activity 
in lattice-QCD studies of the electromagnetic form factors. Significant efforts to probe 
baryon electromagnetic structure in lattice QCD have been driven by the Adelaide g roup 

IF 



13, the Cyprus group 0, and the QCDSF ill and LHP Collaborations Blls 



While lattice-QCD provides the strongest tool for studying non-perturbative phenomena in 
QCD, it does come with its own challenges. In particular, artifacts of (unitary) lattice 
simulations arising from discretisation, finite-volume and unphysical quark masses all need 
to be carefully dealt with to extract the predictions of QCD relevant to the real world. 
(There are further, more complicated issues in dealing with non-unitary approximations: 
partial- and full- quenching, 4th-rooting, mixed-action etc.). 

The principle focus of this manuscript is on the quark-mass dependence of baryon charge 
radii and the chiral extrapolation of lattice simulation results performed at unphysically 
large quark masses. Characterizing the quark-mass dependence of hadronic observables in 
QCD can be achieved within the low-energy effective theory of QCD, xPT — see for example, 
Bernard's recent review on baryon phenomena (23 |. 

A celebrated feature of such effective theories is the model-independence of leading nonan- 
alytic contributions to quark- mass expansions 25|. Such chiral logs are a direct consequence 
of the spontaneously broken chiral symmetry of QCD. Neglecting such behavior in chiral 
extrapolations (of even high-quality lattice calculations) can potentially render results which 
have less resemblance of QCD than some of the models discussed above. While the inclu- 
sion of such logs in chiral extrapolations is necessary to maintain QCD in the extraction 
of physical results 26, [13, 28, [2^, 30], this can be a challenging task because of the poor 
convergence properties of the EFT expansion at moderate quark masses jsil, [H, 33, 3 



In this manuscript we work with finite-range regularisatio n ( FRR) to improve the con- 
vergence properties of the quark-mass expansion of the EFT [3l|, [3^ 
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As an effective theory describing long-distance phenomena, xPT can also describe the 
finite-volume effects of restricting the underlying theory to finite boundary conditions 36 . 
The accuracy of these corrections is also subject to the usual conditions, being at light- 
enough quark masses and large-enough volumes. In this work we impose the finite boundary 
conditions on the relevant one-loop graphs. These provide the leading estimates of the finite- 
volume corrections. 

The lattice simulation results studied in this paper come from the CSSM Lattice Collab- 
oration [17]. These lattice results have been evaluated with quenched gauge-field ensembles. 
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With the underlying dynamics modified by neglecting the quark loops of the QCD vacuum, 
the effective low-energy theory is described by quenched xPT (Q^PT). For baryons, this was 
first formulated by Labrenz and Sharpe [37| and for relevant work on the electromagnetic 



form factors in quenched and partially-quenched theories see Refs. [38|, |39|, |40|, |41 . 

Given quenched lattice results, we use the leading meson-loop diagrams to estimate the 
corrections in going to the fully dynamical theory |32| . This is based on the empirical obser- 
vation that the discrepancies between the quenched and dynamical nucleon and Delta-baryon 
masses are well described by the associated leading meson-loop dressings, as evaluated with 
a dipole finite-range regulator [42]. The physical picture drawn from these results is quite 
intuitive, once the interquark forces are matched at an intermediate distance scale (in this 
case the Sommer scale 43|) the differing long-range features are dominantly described by 
the low-energy EFT. In this case the FRR scale acts to separate the long from the short. 
Further, any residual difference in the qq force at quite short distances does not appear to 
play a significant role in the bulk, low-energy structure. 

In section II, we will briefiy introduce the chiral Lagrangian which is used for the octet 
charge form factors. Charge form factors and radii are calculated in section III with quenched 
and full QCD. Numerical results are presented in section IV and finally section V is the 
summary. 



II. CHIRAL LAGRANGIAN 

There are many papers which deal with heavy baryon chiral perturbation theory. For 



details see for example, Refs. |24j . |44J . |45| . For completeness, we briefiy introduce the 
formalism in this section. In heavy baryon chiral perturbation theory, the lowest chiral 
Lagrangian for the baryon-meson interaction which will be used in the calculation of the 
octet-baryon charge form factors is 

= iTiB^iv ■ V)B, + 2DTrB,Si:{A,„ B,} + 2FTrB,Si:[A^, B,] 

-if^{v ■ V)n^ + C{fi:A^B, + B.A^T^), (1) 

where S"^ is the covariant spin-operator defined as 

S'i = ll'o'^'^v^. (2) 

Here, v'^ is the nucleon four velocity (in the rest frame, we have v'^ = (1,0)). We incorporate 
the explicit propagation of octet and decuplet baryon states, with D, F and C denoting 
the relevant meson-baryon couplings (which, in principle, are to be determined in the chiral 
limit). The chiral covariant derivative is written as D^By = d^By + [V^,!?^,]. The 
pseudoscalar meson octet couples to the baryon field through the vector and axial vector 
combinations 

V, = liCd.C^ + C^d.C), A, = ^iCd.C^ - c^d.C), (3) 



where 



e'^/f, / = 93 MeV. (4) 
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The matrix of pseudoscalar fields (f) is expressed as 




-Bt, and are velocity-dependent fields which are related to the original baryon octet and 
decuplet fields B and by 

B^{x) = e'"*^'^"''^''S(x), (6) 

T^(x) = (^^N^^^^'^T^^x). (7) 

In the chiral 5*^7(3) limit, the octet baryons are degenerate. As the physical strange-quark 
mass is significant, we incorporate the physical mass-splittings in the evaluation of the 
chiral-loop diagrams. 

The explicit form of the octet-baryon matrix is written as 

S= S- -^SO + ^A n\. (8) 

Vs- so --|Ay 

The baryon decuplets are defined by the rank-3 symmetric tensor, with unique elements 
given by 

1 1 

^111 = ^'^'^ ^ Tii2 — -^122 = ' -^222 = ^ , 

t„3^-Le-, r..3^-LE..», r.3^-LE..-, 

ri33 = -^H-'°, r233 = -LH-'-, r333 = S2-. (9) 

In the heavy-baryon formalism, the propagators of the octet and decuplet baryon, j, are 
respectively expressed as 

and , (10) 



vk- A^^ + ie v-k- A^^ + ie ' 

with Pf^" = v'^v" — g'^" — {4/3)3^3'^. A"* = rrih — rua is the mass difference of between the 
two baryons. The propagator of meson k {k = ii, K, rj) is the usual free propagator, i.e.: 

(11) 



k"^ — ml + is 



III. CHARGE FORM FACTORS AND RADII 

In the heavy-baryon formalism, the baryon form factors are defined as 

< B{p')\J,\B{p) u{p') \v,Ge{Q') + '-^^^^^^^^^^^Gm{Q')]u{p), (12) 
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where q = p' — p and = — g^. In this paper, we focus on the charge form factors since we 
will extrapolate octet-baryon charge radii. 

With the Lagrangian, the diagrams for the charge form factors are shown in Fig. [H Our 
counting is ordered in terms of powers of the quark mass contributing to the given electric 
radii. For the purposes of counting, where the typical baryon mass splitting is small relative 
to the lattice meson masses, we treat the power of the quark mass that defines the order of a 
given diagram by neglecting the mass splitting, A. The diagrams are nevertheless evaluated 
with the splittings restored. In the limit << A, diagrams a with the intermediate nucleon 
state and c give rise to the leading-log divergence of charge radii in full QCD. Diagram b 
gives a next to leading order nonanalytic term associated with the mass difference between 
octet and decuplet baryons. For the intermediate hyperson states of diagram a, the mass 
difference between hyperons and nucleon also provides a next to leading order nonanalytic 
term as diagram b. In this limit, we only have a partial contribution of next order. This 
is certainly one of the places where our analysis is somewhat phenomenological. In the 
quenched case, each diagram has a different coefficient from that of dynamical QCD. In 
particular, diagram c has no contribution in the quenched case. Further, the double hairpin 
diagram d contributes only in the quenched case, where the t] is degenerate with the pion. 




1 C 



Octet 

Decuplet 
meson 

y 

FIG. 1: Leading and next to leading order diagrams for the baryon charge form factors. Diagram 
d contributes only in the quenched case. 



Upon integration over k^, the contribution to charge form factors of Fig. la can be 
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expressed as 



M(r^2^_ f^t u{k)u{k - q)k ■ {k - q) 



16^3^2 y uj{k)uj{k - q){uj{k) + uj{k - q)) 



Here, uj{k) = \/ m? + P is the energy of the meson. We note that a summation over 
the relevant intermediate states is assumed. In our calculation we use the finite range 
regular isat ion and is the ultra-violet regulator. Both pion and kaon loops are included 
in the calculation. In the kaon case, Eq. (fT3|l must of course, include the mass differences in 
the intermediate states. In full QCD, the coefficients are determined from the Lagrangian. In 



the quenched case, the coefficients are obtained as in Ref. j39| using the quark flows of Figs. [2] 
and [HI The results are the same as those extracted within the graded symmetry formalism. 
In Fig. [31 the diagrams with 7r° loop are not shown since they have no contribution in the 
quenched or full QCD. They do contribute to the full QCD valence sector and are included 
in our calculation for the valence sector results. 
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The contribution to the charge form factors of Fig. [TJd can be written as 



16^3/2 



dk- 



u{k)u{k — q)k ■ {k — q) 



{uj{k) + A)(a;(fc -q) + ^){uj{k) + uj{k - q)) 



(14) 



where A is the positive mass difference between octet and decuplet baryons. The contribu- 
tion to the form factors of Fig. Ic is expressed as 



167r3/2 



dk- 



u{ky 



uj{k + q/2) + uo{k - q/2) 



(15) 



In the above equations, depends on the baryon type (or quark type), meson loop type 
and quenched or full QCD in the calculation. 

In the quenched case, the double hairpin term from the rj' is expressed as 



967r3/2 



dk^^ ^ Ge{Q')G%^ 



(16) 



where Mq is the double hairpin interaction strength. As a vertex renormalisation in the 
heavy-baryon limit, the Q^-dependence f actor izes to define a Q^-independent G%. We should 
mention that at the lowest order, the double hairpin diagram is independent. The higher 
order terms arising from the dependence of the contributions of this graph at the masses 
probed in the lattice simulations are negligible. 



TABLE I: Coefficient of proportionality for each diagram of Fig. [2] for the proton. 



(a) 



(b) 



(c)(e)(i) (d)(f)(h) 



(g) 



2F' 



-{D-Ff \D^ + 2F^ {D-Ff \{D + Ff 



5^)2 



+ 2DF 



TABLE IL Coefficient of proportionality for each diagram of Fig. [3] for the proton. 



(b) 



(d)(h) 



(e)(i) 



(f)(j) 



TABLE III: Coefficients for quarks in the octet-baryons in full and quenched QCD for Fig. la. 
The intermediate meson is vr. 
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FIG. 3: Feymann diagrams of Fig. Ic for the proton, in terms of quark lines. 



In the above formulas, the coefficients in quenched and full QCD can be obtained following 
the methodology of Ref. 3^. For example, the diagram Fig. la is shown in detail with quark 



lines in Fig. [21 For the pion loop, in full QCD Fig. 2a and Fig. 2c make contributions, while 
in the quenched case Fig. 2a and Fig. 2b make contributions. Therefore, we need to get the 
coefficient for each diagram. The coefficient for Fig. 2c is the same as Fig. 2e which is known 
from the Lagrangian, since QCD is flavor blind. By subtracting this known coefficient from 
the total coefficient of full QCD, we can get the quenched one of Fig. 2a. As an example, 
the coefficient of each diagram of Fig. [2] is listed in Table I. In the same way, the diagram 
of Fig. Ic can be shown in detail in Fig. [3l Table II gives the coefficient of each diagram of 
Fig. [3] for the proton. We can see that the coefficients for the first two diagrams Fig. 3a and 
Fig. 3b are zero which means Fig. Ic has no contribution to the proton charge form factor in 
the quenched case. In fact, for other octet baryons, this diagram has no contribution either 
in the quenched case. 

One can also concentrate on each quark contribution to the form factors. All the coeffi- 
cients for each quark are shown in Tables III- VIII. In these Tables, only three baryons are 
listed. The coefficients for neutron, S~, S~ and A can be obtained by the following 
charge symmetry relations: 
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Ma = 



1 

" 3 



Ms- = (i=o, (is- = Uso, Ss- = -5=0; 

^(m„ + Mho), dso = ^(rfp + d^-), s^o = 



(18) 
(19) 

(20) 



Up + dp + Mho + d^o — 2 + c^s-) 



SA = -(2Sp + 2sH0 -ss+). (21) 



We express our charge form factors as 

1 



(22) 



where Z is the wave function renormahzation constant of charge expressed as Z = GEiQ"^ 



0) ~ YTi=a^^E iQ'^ — 0)- Ge{Q'^ = 0) is the charge of the baryon. G^'' is expressed in 
Eqs. f fT3|) - f|T6|) . Therefore, GEiQ"^) can be written as 

d 

Ge{Q^) = {Ge{Q^ = 0) - 5^Gg(g2 = 0) - -(ao + a^ml + a,ml)Q^ + G^^\Q^) 



(i) 



+Gf(Q^) + G^^^(g^))/(l-G^ 



ic)i 



(23) 



In anticipation of considering charge radii, we have defined the contribution from the double- 
hairpin. Fig. Id, to be proportional to the renormalised form factor using the factorization 
defined in Eq. (|T6|) . This follows a similar procedure to that outlined in Ref. 32 . 



The expansion in characterizes the non-chiral quark-mass dependence of the electric 
charge radius of each baryon. To the leading order we work in this manuscript, ao acts as 
a counter-term to the loops and thereby removing scale-dependence of the formal expan- 
sion. Previous works have shown that including a partial contribution from the next higher 
analytic order beyond which one is working (e.g. Refs. [3l|, H^) can reduce the dependence 
on the regularisation, motivating the a2 parameter. We have also included 04, which has 
simply been included to better describe the lattice data, such that its justification is purely 
empirical (and not mathematical). Its presence mirrors the success obtained with a similar 
approach for the octet baryon magnetic moments. Thereby our form should only be seen to 
incorporate the leading-logarithmic behavior of the EFT with the associated counterterm. 

The Sachs charge radius is defined by 



< >e= -6 



dGEiQ' 
dQ^ 



(24) 



From the expression of Ge^Q"^), < r"^ >e can be written as 
< >E= (o-o + a2ml + a4m^ - 6 



dQ^ 



|q.=o)/(1-Q). (25) 



The above free parameters ao, a2 and a4 are to be determined by fitting quenched lattice 
results with the described quenched loop integrals. The octet charge radii was investigated 
in Ref. [H, 47 1 where the SU{3) symmetry was applied. Our approach is based on an 



SU(2) framework where the strange quark mass is held fixed, and the light quarks are 
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always degenerate. Thereby each baryon isospin multiplet will carry independent LECs. 
Effectively there will be two (independent) Oj for each multiplet. We do not impose the 
symmetry breaking patterns of SU{3), as we regard the strange quark mass as a scale that 
cannot be described well by this low-order expansion. We note it could be interesting to 
compare our effective SU (2) LECs with those extracted from earlier studies. We do not do 
this, but we have now included a comparison of the predictions for the radii of the physical 



states. In Ref. j48|, the authors have developed two-flavor xPT to describe hyperons which 



are embedded into SU{2) multiplets. 



IV. NUMERICAL RESULTS 

In the numerical calculations, the parameters are chosen as D = 0.76 and F = 0.50 
{gA = D + F = 1.26). The coupling constant C is chosen to be —2D, as estimated by SU{6) 
relations — which gives a similar value to that obtained from the hadronic decay width of 
the A. 

Here the finite-range regulator is chosen to take the dipole form 

= (1 + P/AT ' 

with A = 0.8 ± 0.1 GeV. This selected range of A for the dipole has been found to give good 
quantitative estimates of quenching artifacts for baryon masses jl^l and magnetic moments 

Using a fixed strange-quark mass, we estimate the i^-meson mass to obey the relationship 

2 ^ 2 , 2 I ■'"2 1 frt'7\ 



We first study the u-quark contribution to the proton charge radius. Four kinds of 
extrapolations are shown in Fig. HI The square, rhombus, triangle, and round symbols are 
for the finite volume quenched QCD, infinite volume quenched QCD, valence sector and 
full QCD results, respectively. We remind the reader that the "valence" result denotes the 
connected current insertions in Full QCD. The quenched lattice results are described very 
well. For infinite volume, all of them have log-divergent behavior at m,r = 0, which means 
that the pion cloud extends to infinity for a massless pion. For finite volume, the integration 
is replaced by the summation of the momentum which shows no log divergence in the chiral 
limit. Technically, one should be cautious when considering the finite volume curve in the 
domain where m^rL becomes small. Ideally, one would like to stay in the regime of m.,^L > 6 
(or 27r). Ambitious lattice calculations push this down to 4 or 3, at which point one may 
be approaching the limits of describing the finite-volume corrections by the one-loop EFT. 
This limit, for our 2.56 fm box, is at pion masses of order 250 to 300 MeV, the lightest pion 
mass considered in the lattice simulations. This is not of serious concern, as corrections are 
performed at finite volume and extrapolations subsequently performed at infinite volume. 

Fig. [5] displays the d-quark contribution to the proton charge radius. The lines with 
different types have the same meaning as in Fig. 4. For infinite volume, the charge radius 
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FIG. 4: The contribution of a single u quark with unit charge to the proton charge radius versus 
pion mass. The square, rhombus, triangle, and round symbols are for the finite volume quenched 
QCD, infinite volume quenched QCD, valence sector and full QCD results, respectively. 
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FIG. 5: The contribution of a d quark with unit charge to the proton charge radius versus pion 
mass. The square, rhombus, triangle, and round symbols are for the finite volume quenched QCD, 
infinite volume quenched QCD, valence sector and full QCD results, respectively. 



bends down as approaches zero, indicating that the d contributes more strongly to the 
long-ranging tail than the d in accord with the predominant p — *• nvr channel. At m,r = 0, 
the curves also have log divergent behavior. We note that the positive charge radius is 
much like a constituent- quark type expectation, and that the unusual chiral features are 
only anticipated to become dominant far below the physical quark mass. 
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FIG. 6: The contribution of a single u quark with unit charge to the charge radius versus pion 
mass. The square, rhombus, triangle, and round symbols are for the finite volume quenched QCD, 
infinite volume quenched QCD, valence sector and full QCD results, respectively. 
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FIG. 7: The contribution of an s quark with unit charge to the E"*" charge radius versus pion 
mass. The square, rhombus, triangle, and round symbols are for the finite volume quenched QCD, 
infinite volume quenched QCD, valence sector and full QCD results, respectively. 



We now discuss the quark distribution radius in strange octet baryons. The w-quark 
contribution to the S"*" charge radius is shown in Fig. [61 In contrast to the proton case, in 
quenched QCD there is no pion loop contribution to the S"*" charge radius. The non trivial 
loop contributions involve the K meson. Consequently, the w-quark distribution radius has 
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TABLE IV: Coefficients 13% for quarks in the octet-baryons in full and quenched QCD for Fig. la. The intermediate meson is K. 
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TABLE V: Coefficients for quarks in the octet-baryons in full and quenched QCD for Fig. lb. 
The intermediate meson is vr. 
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TABLE VL Coefficients for quarks in the octet-baryons in full and quenched QCD for Fig. lb. 
The intermediate meson is K. 
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no log divergence (as the SU{2) chiral limit is approached). In full QCD, both the valence 
and total sectors exhibit log-divergent charge radii. At the physical pion mass, the total 
M-quark distribution radius is found to be very similar to that in the proton. 

The singly-represented quark of the S"*" is a strange quark. We show its contribution to the 
S+ charge radius in Fig. [71 Since there cannot be any leading order pion-loop contributions 
to the s-quark radius the chiral corrections are much less dramatic than in the light-quark 
sector. One may also compare K contributions at the physical point with those of the light- 
quark sector as observed near ~ "^l'{phys) ~ ^.25 GeV^. Because the s quark mass is 
held fixed, any variation is due to an environment effect associated with the light quarks. 



The is composed of two strange quarks and one up quark. The strange quark contri- 
bution to 'EP is shown in Fig. [HI Similar to the s quark in the S"*", all the curves display a 
very mild environment dependence on the light-quark mass. 

The u-quark contribution to the S° is shown in Fig. [9] In full QCD, both the full and 
valence contributions are divergent as — ^ 0. Here we see that the charge radius diverges 
in the positive direction, opposite to that of the d quark in the proton. Here, the pion loop 
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TABLE VII: Coefficients for quarks in tlie octet-baryons in full and quenched QCD for Fig. Ic. 
The intermediate meson is vr. 
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TABLE VIII: Coefficients for quarks in the octet-baryons in full and quenched QCD for Fig. Ic. 
The intermediate meson is K. 
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3 



must contain a light anti-quark from the sea and thereby the valence u is pure quark (that 
is, it cannot be an anti-quark, or it cannot couple to a pion through a Z-type diagram). 

The single light quark and strange quark contribution to the A is shown in Fig. [10] and 
[TTl respectively. It is interesting that the light quark sector does not show much curvature. 
This is consistent with the chiral coefficients for the leading-order pion dressing vanishing. 
The contribution from the vr"*" and vr" meson clouds cancel each other. The s quark in the 
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FIG. 8: The contribution of a single s quark with unit charge to the charge radius versus pion 
mass. The square, rhombus, triangle, and round symbols are for the finite volume quenched QCD, 
infinite volume quenched QCD, valence sector and full QCD results, respectively. 
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FIG. 9: The contribution of a u quark with unit charge to the H'^ charge radius versus pion mass. 
The square, rhombus, triangle, and round symbols are for the finite volume quenched QCD, infinite 
volume quenched QCD, valence sector and full QCD results, respectively. 
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FIG. 10: The contribution of a single light quark with unit charge to the A charge radius versus 
pion mass. The square, rhombus, triangle, and round symbols are for the finite volume quenched 
QCD, infinite volume quenched QCD, valence sector and full QCD results, respectively. 



A couples strongly to NK and it looks like there is some environment dependence. 

Using the charge-symmetry relations above, Eqs. ( fT711 -( |2T]) . one can reconstruct any 
desired baryon form factors by applying the appropriate charge factors. 

We now discuss the constructed octet-baryon charge radii. The proton charge radii are 
illustrated in Fig. [121 At the physical pion mass, we find the proton charge radius is about 
0.69 ±0.05fm^. The m, d and s quarks all provide a contribution to the proton charge radii. 
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FIG. 11: The contribution of an s quark with unit charge to the A charge radius versus pion mass. 
The square, rhombus, triangle, and round symbols are for the finite volume quenched QCD, infinite 
volume quenched QCD, valence sector and full QCD results, respectively. 
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FIG. 12: The proton charge radius versus pion mass. The square, rhombus, triangle, and round 
symbols are for the finite volume quenched QCD, infinite volume quenched QCD, valence sector 
and full QCD results, respectively. 



We note that the valence sector alone contributes a large fraction to the total charge radius. 
There is a small correction by adding all the disconnected contributions of the 3 light-quark 
flavors in the process of correcting QQCD via FRR EFT. Further, the strangeness component 
only constitutes a small part of this small correction HO • 



The charge radius of S"*" is shown in Fig. [131 The numerical value is observed to be a little 



larger than that of proton, with the physical charge radius of the S"*" being 0.75 ± 0.05 fm 
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The primary cause of this radius being larger than the proton is the fact that the (negatively- 
charged) strange quark in the S has a narrower spatial distribution than the down quark in 
the proton. 
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FIG. 13: charge radius versus pion mass. The square, rhombus, triangle, and round symbols 
are for the finite volume quenched QCD, infinite volume quenched QCD, valence sector and full 
QCD results, respectively. 




In Fig. [HI the charge radii of the octet baryons at the physical pion mass are shown. 
The extrapolated physical radii of the proton, neutron and E~ are in good agreement with 
the experimental data. The radius of H° is positive, in contrast to the neutron radius. The 
reason is that, for the neutron, in full QCD, the vr" cloud is the dominant contribution, 
whereas for the the pion intermediate state is a 7r+ (with a S~ intermediate baryon). 
Since both tt"*" and vr" contribute to the A, the radius of the A is close to zero. This is perhaps 
surprising as one would expect the net positive charge of the light quarks to dominate more. 

The results of our extrapolation and unquenching are summarized in Table IIXI Here, 
we see a careful breakdown of the propagation of uncertainties at the various stages of 
calculation. 



V. SUMMARY 



We extrapolated state-of-the-art lattice results for the quark-sector decomposition of 
the octet-baryon charge radii in quenched heavy baryon chiral perturbation theory using 
FRR. All leading-loop diagrams have been incorporated, including all contributions from 
both octet and decuplet baryon intermediate states. Finite-range regularisation has been 
utilized in the one loop calculation to improve the convergence at moderate quark masses. 
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FIG. 14: Octet-baryon charge radii at the physical pion mass. The hollow square, solid rhombus, 
triangle and square symbols are for the finite volume quenched QCD, infinite volume quenched 
QCD, valence sector and full QCD results, respectively. The experimental data for proton, neutron 
and S~ is shown with the left-most bullet. 



TABLE IX: Summary of extrapolation and unquenching results. Each major column consists of 
our best value (in fm^), followed by the three dominant sources of uncertainty: statistical, lattice 
scale determination a = 0.128 it 0.006 fm and regularization scale A = 0.8 ib 0.1 GeV, respectively, 
all quoted relative to the final digit of the best value. The sign of the uncertainty reflects the 
correlation with a or A. 
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0.584 


46 


45 


16 
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Further, the use of FRR provides a separation of scales, which has enabled the use of 
a demonstrated technique to obtain estimates of the full-QCD results from the quenched 
lattice simulations. We acknowledge the phenomenological aspects of our calculation. As 
such, we have an empirically motivated extrapolation form which is equivalent to the leading 
one-loop HBChPT form in the domain where higher-order terms are genuinely negligible. 

The individual quark contributions to the baryon charge radii in quenched QCD have 
been extrapolated, with subsequent predictions for the corresponding radii in both valence 
and full QCD. We note that the valence predictions will be readily confronted with the 
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next generation of lattice simulations of full QCD, as only connected current insertions 
are required. The contribution from pion loops is observed to generate significantly more 
enhancement of the radii than that from K mesons. Further, disconnected contributions to 
the radii are rather small compared with the dominant valence contributions. One of the 
consequences of this is that the strange-quark disconnected contributions only represent a 
small component of an already small correction. 

The charge radii of the proton, neutron and are in good agreement with the exper- 
imental results. Our result for E~ is quite accurate relative to the current experimental 
measurement. Further, we have presented predictions for the remaining five baryons, for 
which no experimental information exists to date. These predictions remain to be tested by 
potential future measurements and future studies in ab initio lattice QCD. 
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